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Abstract. 

We present high-resolution echelle spectroscopy, obtained with the UVES spectrograph on ESO/VLT, of two 
luminous star clusters in the metal-poor blue compact galaxy ESO 338-IG04 at a distance of 37.5 Mpc. Cross- 
correlating with template stars, we obtain line-of-sight velocity dispersions of 33 and 17 km s _1 . By combining 
with size estimates from Hubble Space Telescope images we infer dynamical masses of 1.3 x 10 7 M Q and 4.0 x 10 6 
Mofor the two clusters, making them among the most massive known. The less massive cluster is the faintest 
cluster for which a dynamical mass has yet been obtained. In both clusters we detect Balmer absorption lines 
which we use to estimate their ages. From the younger (~ 6 Myr) and more massive cluster, we detect Hell 
A4686 emission of intermediate width, indicating the presence of very massive O-stars. Moreover, analysis of the 
[O ill] A5007 and Ha emission lines from the region near the younger cluster indicates that it is associated with a 
bubble expanding at ~ 40 km s _1 . We also see from the Na iD absorption lines indications of neutral gas flows 
towards the younger cluster. We compare the dynamical masses with those derived from photometry and discuss 
implications for the stellar initial mass function (IMF) in each cluster. Both clusters are compatible with rather 
normal IMFs which will favour their long-term survival and evolution into massive bona fide globular clusters. 

Key words, galaxies: evolution - galaxies: individual: ESO338-IG04 (= Tol 1924-416) - galaxies: starburst - 
galaxies: star clusters - galaxies: stellar content 

1. Introduction such photometric masses have been determined, only a 

small fraction of extragalactic SSCs and GCs have dynam- 

Super star clusters (SSCs) are important sites fo r star ^ masg estimates . Taken together, measurements of the 

formation in starburst galaxies (Arp and Sandage photometric and dynamical masses can allow constraints 

Meurer et al. 1995), and may be the progenitors of globu- Qn the gtdlar initial magg function (IMF) and the dynam . 

lar clusters (GCs) . Different observational techniques have ical gtate in thege objects 

been used to investigate the nature and history of these The spectrum of a cluster can provide other informa . 

objects. High resolution imaging, notably with the Hubble ^ aW ^ gteUar population . For example, Gonzalez 

Space Telescope (HST), has been used to collect photom- Delgado et &l inferred & population of Wolf-Rayet 

etry and to estimate the sizes of many SSCs (Whitmore starg ^ the masgive clugter NQC 15m _ A ^ and 0riglia ^ 

2003). Photometry, m combination with spectral evolu- ^ j^j uged u i travio i et spectra to show that the same 

tionary synthesis models, has also beenused to estimate dugter hag a population of - S tars. 

cluster ages and masses (Ostlin et al. [Mi . While many High-resolution spectra can be used to measure the ve- 

.. locity dispersion, and, assuming virialisation, to estimate 

Send offprint requests to: G. Ostlin j ' • i c oc/~i j /"in tt j m-i- i 

* , , . _ „ dynamical masses for SSCs and GCs. Ho and Mlippenko 

Based on observations collected at European Southern ., ,,,,,, ■ ■ _, nn „t ■ „ , , , . 

„, Ti i -i j u ' rr ivt 1996a, 1996b were hrst to use the cross-correlation 

Observatory, Paranal, Chile, under observing programme 65. JN- ' ' 

0704; and on observations with the NASA/ESA Hubble Space ^od described by Simkin and Tonry and Davis 

Telescope, obtained at the Space Telescope Science Institute, <X2Z9 to measure the velocity dispersions of NGC 1569- 

which is operated by the Association of Universities for A (giving a kinematic mass of 3.3 x 10 5 M Q ) and the 

Research in Astronomy, Inc., under NASA contract NAS5- SSC in NGC 1705 (8.2 x 10 4 M ). The same meth od has 

26555. been applied to M82-F by Smith and Gallagher (2001 
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Fig. 1. Acquisition. Upper left: HST/WFPC2 image of ESO 338-IG04 in the F814W band, rotated to match the po- 
sition angle used for our UVES observations. The orientation is indicated and the two main targets for the spectroscopy 
discussed in this paper are labelled. Upper right: Enlargement of the upper left image showing the slit location, the 
main targets and some other nearby clusters. The width of the slit is 0.7". The clusters have been identified with the 

mm 



index number used by Ostlin et al. (1998). All but #34 belong to their 'inner sample'. Lower left: The image in the 
upper left panel, convolved with a Gaussian kernel the same size as the seeing (FWHM=0".65) to simulate the actual 
observing conditions. Lower right: Slit view image from the red camera of UVES. The slit is the horizontal bright 
bar, and has a length of 11". The bright foreground star on the right was used for offsetting in the acquisition. 



1.2 x 10 6 M ), to five clusters in NGC 4214 and NGC 
6946 by Larsen et al. (jMUl IM3I 0.2 - 1.8 x 10 6 M ), 
by Gilbert and Graham (|2002[1 to clusters Al, A2 and B 
in NGC 1569 (3.9 x 10 5 M Q , 4.4 x 10 5 M Q and 2.3 x 10 5 



(2003) to two SSCs in M 82 (1.5 x 10 6 M Q and 3.5 x 10 5 
M Q ). Mengel et al. l|2"0"02) measured the dynamical masses 
of five SSCs in NGC 4038/4039 (0.65-4.7x 10 6 M Q ) by fit- 
ting model spectra in the near-IR and optical. Maraston et 



M Q , respectively) and by McCrady, Gilbert and Graham al. (2004) presented a dynamical mass estimate of as much 
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as 8 x 1O 7 M for the cluster W3 in NGC 7252, which may 
possibly be a dwarf galaxy. Bastian et al. (2006) found 
masses of 1.6 x 10 7 M Q for both W30 in NGC 7252 and 
G114 in NGC 1316, and present a compilation of previ- 
ously published results (their Table 6). 

A possibly related class of objects are the nuclear 
star clusters found in many bulgeless spiral galaxies (e.g. 
Walcher et al. I2005j) . some of which appear to be very 
massive but which probably have more extended star for- 
mation histories than GCs and SSCs. 

ESO338-IG04, also known as Tol 1924-416, is a lumi- 
nous (Mb = —19) blue compact galaxy (BCG) — the 
closest in a class of galaxies that are rare in the local 
universe but that become increasingly important towards 
higher rcdshifts and the peak in the cosmic star forma- 
tion rate at z ~ 1 (e.g. Werk et al. l2lTO|l . ESO 338- 
IG04 shows both vigorous star formation and a rich pop- 
ulation of more than 100 young SSCs and old globular 
clusters (Ostlin et al. I1998I1 1 . Analysis of the age dis- 
tribution of these clusters shows evidence that at least 
one strong starburst occurred in the past, a couple of 
Gyr ago, and the present starburst has probably been ac- 
tive for about 40 Myr (Ostlin et al. lifoTj3"|l . In addition 
ESO 338-IG04 hosts a population of cosmologically old 
(age < 10 Gyr) globular clusters. What makes ESO 338- 
IG04 particularly interesting as a star cluster formation 
laboratory is its combination of small reddening, many 
clusters, and a low metallicity (12 + log(0/H) = 8.0). 
Compared to other metal-poor star-forming galaxies in the 
local universe (such as NGC 1569, NGC 1705, and He 2- 
10) the clusters in ESO 338-IG04 are more luminous and 
numerous. The present starburst has been triggered by 
a small merger or interaction with its companion galaxy 
(Ostlin et al. l2DCjTl Cannon et al. 2004). Hence ESO 338- 
IG04 offers the possibility of directly studying globular 
cluster formation associated with hierarchical galaxy evo- 
lution, at a low metallicity. 

In this paper, we report high-dispersion echelle spec- 
troscopy of two of these clusters, labelled Inner-#23 and 
Outer-#34 by Ostlin et al. I|1998(l . Hereafter, we refer to 
them simply as #23 and #34. Cluster #23, a young blue 
cluster situated in the centre of the starburst region, is 
the most luminous in the whole galaxy (My = —15.5). 
Cluster #34 is of intermediate age, ~ 1 Gyr, but still very 
luminous (My = -12.8), which led Ostlin et al. (|E)9"8"|) to 
suspect a mass in excess of 10 7 Mq. Here, we use high- 
dispersion VLT/UVES spectroscopy to make dynamical 
mass estimates and new age estimates for both clusters. 
These two clusters were selected since they represent the 
most luminous of the young and intermediate age clusters, 
respectively, and could be observed in a single UVES slit. 

The outline of the rest of the paper is as follows. In 
Sect. 2 we describe the observations and reductions. In 
Sect. 3 we show the results obtained; subsection 3.4 in 
particular describes how we obtained estimates of the ve- 

1 Following Ostlin et al. 119988 we adopt a distance of 37.5 
Mpc to ESO 338-IG04. 



locity dispersions of the clusters through cross-correlation 
techniques. In Sect. 4 we discuss the inferred dynamical 
masses and compare with results derived from photome- 
try. Sect. 5 summarises our conclusions. 

2. Observations and reductions 

2.1. UVES-spectroscopy 

Our spectra were taken on 2000 June 15 using the echelle 
spectrograph UVES mounted on the telescope Kueyen 
(UT2) of the Very Large Telescope, at the European 
Southern Observatory's Paranal site. The observations 
were performed in service mode. The seeing was typically 
0".6. 

We placed the slit at a fixed position angle of 42°. 6 
in order to capture the spectra of both cluster #23 and 
#34 in the same exposure (see Fig.^l. To minimise con- 
tamination from other sources in this crowded starburst 
environment, we used a slit width of 0".6. The slit lengths 
used were 8" and 11" for the blue and red cameras, respec- 
tively. The atmospheric dispersion corrector in UVES was 
used to minimise slit losses due to differential refraction. 

Using the dichroic setting DIC1 with cross-dispersers 
CD2 and CD3, we covered the optical spectrum from 
3280 A to 6660 A at a resolution of R ~ 60000. We binned 
the CCD pixels 2 x 2 to lower the influence of the readout 
noise, which resulted in a spatial scale of 0".50 for the blue 
arm and 0".36 per pixel for the red arm. A total of eight 
2000-second exposures were taken. The data were reduced 
using the UVES pipeline, with extraction and sky windows 
modified to allow for the extraction of two non-centred 
sources. Our wavelength calibration, done by comparison 
with ThAr lamps, should be accurate to better than 0.02 
A, or about 1 km s _1 . 

Using the same instrumental set-up, we observed five 
template stars (for details, see Table 0. Two of these, 
HR 7706 and HR 7749, were observed on 2000 April 22; 
the rest were observed on the same night as the cluster 
observations. 

2.2. Source extraction 

Special care was taken with the background subtraction, 
though the small size of the slit and the complex H II re- 
gion emission meant that the subtraction is unreliable in 
the vicinity of strong emission lines. This is also due to the 
intrinsically different spatial distribution of nebular emis- 
sion and the stellar continuum (see Ostlin et al. 120031 and 
Fig. El. Spectra were extracted using average extraction, 
centred on the continuum of the clusters. The extraction 
windows were 3 and 5 pixels wide for the blue and red 
arm spectra, respectively. We extracted the spectra for 
both the clusters and the template stars in this way. 

While clusters #23 and #34 are well-separated on the 
sky, light from a number of other sources was also col- 
lected by our slit (see Fig. In particular, we expect 
a contribution from at least two more clusters from the 
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Table 1. Template stars 



Designation Type and class Radial velocity (km s i ) Reference 1 

HR 6961 K4 lb -2.291 Nidever et al. (p002 1 

HR 7277 Kll +2.8 Ochsenbein RW)) 

HR 7706 Kim +1.0 ±0.3 This work 

HR 7749 F5v -30.8 Kharchenko et al. (2004) 

HR 8299 G5iii -57.1 + 0.4 This work 



Note: (1) All references are to catalogues available at the VizieR website, http://vizier.u-strasbg.fr/ (Ochsenbein, Bauer and 
Marcout ["" 



inner sample of Ostlin et al. (1998): numbers 24 and 27 
(#24 and #27, see Fig. QJ. Both are slightly offset from 
the centre of the slit and we estimate that about 50% 
of the flux from #24 and about 15% of that from #27 
are collected by our extraction window. Both clusters are 
however fainter than #23, by 4.2 and 2.1 magnitudes, re- 
spectively, in the F555W filter (Ostlin et al. lTMSjl . Based 
on this we estimate that on the order of 1% or less of 
the flux from #23 has its origin in each of #24 and #27. 
Contamination from known sources should thus not be a 
problem. 

Nevertheless, since #23 is located in a crowded area, 
we also tried using narrower extraction windows (width 1 
pixel for the blue arm, and 3 pixels for the red arm) with 
the centre offset by 1 pixel to the north-east, where the 
crowding is less severe (see Fig.^l. 

2.3. HST Imaging 

In addition to the UVES spectra, we used images from 
the Hubble Space Telescope obtained with the Planetary 
Camera (PC) aperture of the WFPC2 in the F218W, 
F336W, F439W, F555W and F814W p assba nds. Thes e 
data are described fully in Ostlin et al. I|1998l and I2UU3(I . 
We used these images to estimate the sizes of the two clus- 
ters, and to constrain their age and photometric mass from 
comparison with spectral evolutionary synthesis models. 
The photometry used in the current paper differs slightly 
from that of Ostlin et al. ( 1998) in that we have optimised 
the photometric parameters specifically for these two clus- 
ters and also corrected for the charge transfer (in) efficiency 
of WFPC2 (Whitmore et al. I1999al Dolnhin H-ITO)) using 
the web-based tool CTE Tool #1 2 . This effect is noticeable 
for faint sources on low backgrounds and is here impor- 
tant mainly for the F336W magnitude of #34 where it 
amounts to ^0.06 magnitudes. 



3. Results 

3.1. The spectrum of cluster #23 

Figure E] shows parts of the spectrum of cluster #23. The 
spectrum shows a stellar continuum with superimposed 
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Fig. 2. Part of the spectrum of cluster #23 with line iden- 
tifications. The spectrum has not been flux-calibrated. 
Note the absorptions in the Balmer series and the broad 
emission line in Hen A4686. 



narrow emission lines and a few emission lines of inter- 
mediate width. The narrow emission lines come from ex- 
tended emission close to the cluster. Background subtrac- 
tion across this emission was not always successful, lead- 
ing to the spurious absorption features seen in some of the 
emission lines. 

The narrow line emission presumably originates from a 
combination of H II regions and supernova remnants. We 
have identified lines of H I, He I, [N i], [N n], [O n], [O III], 
[Ne in], [S II], [S in] and [Fe m]. The ratio of [O n] A3729 
to A3726 indicates a emitting electron density of around 
150 cm~ 3 . 

In addition, broader emission is seen in He n A4686 
(FWHM ~ 460 kins" 1 ) and possibly C iv A4658 (blended 
with narrower [Fe in]). The line widths we see here are 
comparable to those seen in Of stars, but narrower than 

mm. 



CTE Tool #1 by Andrew Dolphin 



in Wolf-Rayet stars (cf. Nota et al. 

Stellar absorption lines with broad wings are seen in 
the Balmer series. Narrower absorptions are seen in He I, 
Na I D (both at the redshift of ESO 338-IG04 and at zero 



www.stsci.edu/instruments/wfpc2/Wfpc2_cte/wfpc2_cte_calc.htmledshift from the local Galactic ISM) and Ca II K. A num- 
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ber of other absorption features are apparent towards the 
red parts of the spectrum, in particular at 5168/5172 A 
(Mg I), 5185 A (Mg I, MgH), and also at 6138/42 A and 
6496 A. 

3.2. The spectrum of cluster #34 

Cluster #34 was detected with much poorer signal-to- 
noise ratio than #23. Its spectrum shows Na I D absorp- 
tion lines, the Mg I doublet at 5168/5172 A, and a number 
of other features also seen in the template stars. No broad 
emission lines of the type seen in cluster #23 were de- 
tected, though features with the same equivalent widths 
(EW) could easily be hidden in #34's noisier spectrum. 

3.3. Cluster ages from Balmer lines 

We have estimated the ages of the clusters by fitting the 
model spectra of Gonzalez Delgado and Leitherer (1999) 
to the Balmer and He I line spectrum around US and in 
the region from H8 to Hll. 

First we smoothed our spectra to match the resolution 
of the models. Comparing the models with the smoothed 
data, we located the combination of effective temperature 
and surface gravity which gave the best fit to the observed 
line spectrum. Since the line cores are dominated by emis- 
sion in our data, we matched based on the wings of the 
lines. 

For cluster #23 we get T cS = 3.75 ± 0.5 x 10 4 K (cor- 
responding to late O stars), log g — 4.5—5; and for #34 
we get T off = l.Sto 7 3 x 1q4 k ( late B stars), log g = 
4.5—5. In Figures El and 0] we show some of these fits. 
Adopting log g — 4.5, we estimated the equivalent widths 
of the Hi5 and H8 lines using the data in Tables 4 and 
5 of Gonzalez Delgado and Leitherer fiUM GDL99). We 
obtain, for #23, Wm = ^to.l A and Whs = 2.8±^ A, 
and for #34, W m = 12±| A and W m = 9±l A. 

Finally, we estimated the age of each cluster by com- 
paring these equivalent widths derived from H<5 and H8 
with the values for different instantaneous population ages 
presented by GDL99 (their Table 5, for Z=0.001 and a 
Salpeter lfKJ55| l IMF 3 with mass range 1 to 80 M Q ). We 
obtain logt(yr) = 6.8 ± 0.2 for cluster #23, or 6tt Myr. 

This metallicity is slightly lower than the ISM abun- 
dance of ESO 338-IG04, which could lead us to overesti- 
mate the age slightly. Moreover, one may speculate that 
the high dynamical mass of #23 (see Sect. 3.7) could al- 
low some self-enrichment with the result that this cluster 
would be richer in metals than its environment. On the 
other hand, the use of an upper mass limit of 80 M Q , as 
in the GDL99 model, would lead to an underestimate of 
the age if the true IMF were to extend to 120 M Q . 

For the fainter cluster #34, we cannot set very strin- 
gent limits on the age. One reason is the noisy spec- 
trum, leading to uncertain EW estimates. In addition, 

3 We parameterise the IMF as a power law, dN/dM oc M~ a , 
where a Salpeter IMF has a = 2.35 
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Fig. 3. Normalised spectrum of cluster #23, with su- 
perimposed models of Gonzalez Delgado and Leitherer 
(1999). We show models for logg = 4.5 and effective tem- 
peratures of 30 000 K (thin solid line), 35 000 K (thick 
solid line), 40 000 K (broken line) and 45 000 K (dotted 
line). For ease of comparison, the slope of the model spec- 
tra shortward of 3940A has been adjusted to match the 
data. 

the Balmer line absorption EW of a low-metallicity stellar 
population increases with time up to an age ~ 1 Gyr, af- 
ter which it decreases again (Bica and Alloin 1986.). Based 
on the GDL99 models, the estimated EW indicates an 
age ~ 0.3 Gyr, with a lower limit on the age of 50 Myr. 
Comparing also with the empirical US data from Bica and 
Alloin (1986 ), we find for #34 an allowed age range from 
about 0.5 to 3 Gyr. 

As a consistency check, we compared the effective tem- 
perature of the best-fitting models to the effective tem- 
perature at the main-sequence turnoff in the models of 
Bcrtelli et al. I|1994[) . Stars at the turnoff are expected to 
dominate the blue light from the clusters. This estimate 
gives logt(yr) = 6.8 ±0.25 for cluster #23. For #34 we ob- 
tain logt(yr) = 8.6 ±0.6 for a low metallicity (Z = 0.001) 
and ~ 0.5 dex lower for a high metallicity (Z = 0.02). 

Thus we find an age of ~ 6 Myr for cluster #23, and a 
rather generous age span for #34. These results are consis- 
tent with our photometric modelling (Ostlin et al. 2003), 
which we discuss in Sect. 0] 

3.4. Cross-correlation analysis: radial velocities and 
line-of-sight velocity dispersions 

To estimate the radial velocities and line-of-sight velocity 
dispersions of the two clusters, we cross-correlated their 
spectra with a number of template stars, following the 
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Fig. 4. As Figure|3I but for cluster #34. Here the models 
have effective temperatures of 10 000 K (thin solid line), 
15 000 K (thick solid line) and 20 000 K (broken line). 



method of Simkin l|1974|l . Tonry and Davis (1979 ) and Ho 
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Fig. 5. Part of the spectra of the two clusters, showing 
the region around the Mg I triplet, compared to the same 
region in the five template stars. The spectra have been 
smoothed with a 9-pixel boxcar and shifted to a rest wave- 
length scale for clarity of presentation. 



and Filippenko il996alll"996p|) . For this we used the IRAF 4 
task f xcor. 

3.4.1. Selecting the proper template stars 

For the cross-correlation method to yield trustworthy 
measurements of the velocity dispersion, the absorption 
line spectrum should be dominated by stars with intrin- 
sically narrow lines common to both the target and the 
template star. Cool giants and supergiants are likely to be 
the best templates. 

At the age of cluster #34, we expect the luminosity 
to be dominated by late-type giants, which suggests that 
the optimum template star should be of this type. Cluster 
#23 is much younger, but following the discussion in Ho 
and Filippenko (1996a}, we still expect its absorption line 
spectrum in the visual wavelength region to be dominated 
by red supergiants. Examination of the high resolution 
spectra of local stars presented by the UVES Paranal 
Observatory Project (uvespop; Bagnulo et al. 2005|) con- 
firms that in our wavelength range, only stars of spectral 
type later than A are likely to contribute significantly to 
the observed absorption line spectrum. 

As an independent quantitative check, we used 
Starburst99 (Leitherer et al. 119991 see below) to calcu- 
late the contributions from stars of different spectral types 
to some of the well-known Lick indices (Worthey et al. 
1994). Although these are defined for resolutions much 
lower than ours, the type of stars that dominate a partic- 
ular index should also dominate the individual lines that 
make up the index. The cross-correlation method should 
also be sensitive to these lines. Equally, a spectral type 
that gives a negligible contribution to the indices in the 
visual wavelength region is unlikely to contaminate the 
cross-correlation signal. 

At the Starburst99 website 5 it is possible to run simu- 
lations that output the number of stars for each spectral 
type and luminosity class, as a function of age. We ran a 
few such models (see Table [3J| and used this information, 
together with the luminosity and temperature relations 
from Drilling and Landolt (j2000i and the fitting functions 
of Worthey et al. (1994), to predict the relative contri- 
bution to the Lick indices from stars of different spectral 
types as a function of stellar population age. Model Gl 
has the same tracks and metallicity as the standard model 
used in Ostlin et al. (2003) for estimating stellar popu- 
lations in the clusters in ESO338-IG04. Models PI and 
P2 bracket model Gl in metallicity, and also bracket the 
observed nebular abundance for ESO 338-IG04. Finally, 
model P3 has a lower upper mass limit; in Sect. 14.11 we 
discuss the possibility of an IMF devoid of very massive 

4 IRAF (Image Reduction and Analysis Facility) is dis- 
tributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement 
with the US National Science Foundation. 

5 http://www.stsci.edu/science/starburst99/ 
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Table 2. Starburst99 models used to test which stellar 
types dominate the absorption line features in the spectral 
region 4600-6540 A. 



Model 


Tracks 


Z 


M up 
(Ms) 


PI 


Padova 


0.004 


120 


P2 


Padova 


0.0004 


120 


P3 


Padova 


0.004 


20 


Gl 


Geneva 


0.001 


120 



stars. All models use a lower mass limit of 0.1 M and an 
IMF with slope a = 2.35. 

We found that for all ages, only spectral types later 
than A contribute to the indices in the spectral region 
4600-6540 A. For models PI, P2 and Gl and ages younger 
than 4 Myr, these indices are dominated by main sequence 
stars of class F, G and K. Later, red supergiants (RSGs) 
take over completely, and only after 15 Myr do giants and 
supergiants of spectral type earlier than K contribute. If 
the cluster has an IMF lacking very massive stars, as in 
our model P3, the situation is somewhat different. Here 
the main sequence stars of types F to K dominate longer, 
up to ages of just over 10 Myr, after which RSGs take 
over. This is understandable since less massive stars take 
longer to evolve off the main sequence. Hence, if the high- 
mass IMF is truncated at M up = 20M Q , our Fv tem- 
plate star would be the best choice for cross-correlation. 
However, a low-resolution long-slit spectrum of #23 taken 
with VLT/FORS2 (Cumming et al. 2006, in prep.) shows 
the presence of the Ca n near-IR triplet, with an equiv- 
alent width of several A. This indicates that red giants 
or supergiants must be present (Diaz et al. I1989J1 . The 
Starburst99 models also predict that strong Ca n triplet 
absorption lines develop when RSGs appear. 

In summary, even though #23 is young, we can still 
safely assume that red supergiants dominate the absorp- 
tion features we see, and that such stars should be used 
as template stars for cross-correlation. 

Cluster #34 is somewhere in the range 0.3 to 1.4 Gyr 
old, depending on metallicity (see Sect. 13.31 and !4.2|) . For 
such ages, our analysis using Starburst99 indicates that 
the absorption features in the visual wavelength region 
are dominated by red giants of luminosity classes n and 
III, with a contribution from main sequence stars of types 
F to K. Hence a mix of all the template stars in Tabled 
would seem appropriate. 

A remaining potential problem is that the macroturbu- 
lent broadening in our template stars may differ from the 
macroturbulence characteristics of the cluster stars. This 
difference may be of the order of 1-3 km s _1 for giants and 
supergiants fGrav lT^Tl Gray and Toner ITT)87[1 and some- 
what larger for main sequence stars (Gray 1984). Since gi- 
ants and supergiants are likely to dominate the line widths 
in both clusters, the importance of macroturbulence will 
be small, and will not be considered further. 



Table 3. Spectral regions used for cross-correlation in 
this paper in the rest frame of ESO 338-IG04. 'G' and 'R' 
refer to the short and long wavelength detectors, respec- 
tively, of the UVES red arm, avoiding the strongest emis- 
sion lines and telluric absorption. Weaker emission lines 
remaining have been masked out. The subscripts 'B' and 
'R' correspond to the 'blue' and 'red' parts of each region, 
respectively. 



Rest wavelength range (A) Name 


4606-4952, 5023-5552 


G 


5023-5552 


Gr 


5665-6273, 6320-6544 


R 


5665-6273 


Rb 


6320-6544 


Rr 



3.4.2. Selecting the optimum wavelength region 

We used f xcor to carry out separate cross-correlations for 
a number of different wavelength regions. 

Prior to cross-correlation, we masked out strong emis- 
sion lines, telluric absorption features, and absorption in 
the Na I D lines (see discussion in section l3~51 below') . The 
spectrum of #23 shortward of 4745 A is contaminated by 
nebular emission and early-type stars. For cluster #34, 
this part of the spectrum has in addition poor signal-to- 
noise ratio. We therefore omitted this spectral region from 
the cross-correlation analysis. 

We show the names and ranges of the spectral regions 
we have used in Table Inspection of the spectra (see for 
example Figure 0) showed that those features which are 
visible above the noise in the cluster spectra tend also to 
be present in all our template stars, though with different 
strengths. All the template stars gave a positive cross- 
correlation signal between 4606 A and 6545 A, and always 
at the same radial velocity, to within ±4 km s _1 . The blue 
part of region G (4606-4952 A) gave less reliable results, 
and will not be further considered. We found that the 
best cross-correlations were obtained if we filtered out the 
lowest frequencies in the spectrum. For this we used a 
ramp filter rising from zero to unity at wavenumber 90 
(40) for spectral regions G and Gr (Rb, R and R R ). 

For both clusters, the highest cross-correlation signal 
was obtained for spectral region Gr, irrespective of which 
template star was used. Figures and show examples of 
cross-correlation functions (CCFs) for one of the template 
stars, HR 6961. These figures illustrate how the CCF peak 
amplitudes differ between wavelength regions. The CCF 
peak amplitudes for region Gr are always significantly 
larger than regions G, Rb and Rr, which have peak values 
similar to that for R. 

The arguments for the choice of template star pre- 
sented in Sect. 13.4.11 based on the analysis using 
Starburst99, are of course best founded for wavelength 
regions that include some of the Lick indices (Worthey 
et al. I1994[> , This is fulfilled for region Gr, about half of 
whose wavelength range is covered by these indices. The 
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Fig. 6. Examples of cross-correlation functions of cluster 
#23, here using the K4 supergiant HR 6961 as template 
star, in spectral regions R (top) and Gr (middle; see Table 
for details) and in the region around the Na I D lines 
towards #23 (bottom) . The ranges in rest wavelength used 
for each cross-correlation are shown by the horizontal bar 
in each plot: the bar stretches from 4500 A to 7000 A. 



o 



0.2 







0.2 - 







Cluster #34 

i — i — i — i — r 



"i — i — r 



F 



4500 A(A) 7000 





2600 2800 3000 

V hel 10 ( km S_1 ) 

Fig. 7. Same as Figure but for cluster #34 using the 
same template star. The upper panel shows spectral region 
R, the lower panel Gr. 

Table 4. Notation for cross-correlation analysis. 



Symbol 


Definition 




Vheiio 


Heliocentric radial velocity 




w v 


Measured FWHM of the CCF peak 




AW 


Estimated uncertainty of W v 




h 


Peak value of the CCF 




Ah 


Tolerance in h matching 




a 


Velocity dispersion 




Sw 


RMS scatter in W v of derived W v — 


a relation 




RMS scatter in a of derived W v — a 


relation 



other wavelength regions defined in Tableware clearly less 
suitable from this point of view. While region Rb contains 
five defined indices, two are TiO and one is Na D, which 
we have masked out to avoid contamination by interstel- 
lar absorption. The TiO lines could be problematic since 
they are strong only in very cool stars with later spec- 
tral type than any of our template stars, but which may 
nevertheless be present in our clusters (e.g. M giants and 
supergiants). In summary, region Gr appears to be the 
safest to use. It is also the wavelength region which con- 
tains most of apparent absorption features and for which 
the strength of the cross-correlation data is clearly the 
strongest. 



3.4.3. Calibrations and uncertainties 

For each cluster, spectral region and template star, we 
measured the heliocentric radial velocity (fhciio), the am- 
plitude (h) and full width at half maximum (W v ) of the 
best Gaussian fit to the peak of the cross-correlation func- 
tion. To estimate AW, the uncertainty in W Vl we made 
several measurements with fxcor of the width of each 
CCF peak, varying within plausible limits the baseline 
and velocity region for the fit, and taking the standard 
deviation of several measurements. Our notation is sum- 
marised in Tabled 

We measured Uhoiio independently for each cluster and 
template star, using accurate radial velocities from the lit- 
erature (Tabled . Since the catalogue velocity of HR 8299 
was only accurate to 10 km s _1 , and no catalogue veloc- 
ity was available for HR 7706, we cross-correlated these 
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Wavelength region Gr 5023 — 5552 A 
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Fig. 8. An example plot showing how we have estimated 
velocity dispersion a by measuring the width of fits to 
the CCF peak. Diamonds represent models which match, 
to within ±0.015, the observed CCF peak amplitude for 
cluster #23. Crosses match the peak amplitude for #34. 

with HR 6961 and measured new heliocentric velocities 
for these stars (Table QJ. 

To estimate the line-of-sight velocity dispersion a from 
width of the CCF, we followed the method of Ho and 
Filippenko (1996a, 1996b). The basic procedure was as 
follows: we first broadened each template star spectrum 
with Gaussians of known tr, which we then cross-correlated 
with the corresponding observed template star spectrum. 
From the resulting CCF, we measured the width (W v ) 
and amplitude (h) of the best Gaussian fit to the peak, in 
the same way as for the cluster spectra, and used this to 
calibrate a relation that gave the velocity dispersion a as 
a function of the measured CCF width W v . 

Our clusters are however faint, and our CCFs peak at 
lower values than those of previous authors. We therefore 
adapted the method in order to give a more careful cali- 
bration for our data. 

To explore how the quality of the spectra affects the 
W v — a calibration, we added different amounts of ran- 
dom, Gaussian noise to the template star spectra using 
the IRAF task mknoise. Each of these noisy model spectra 
was then cross-correlated with the corresponding observed 
template star spectrum. From among the resulting CCFs, 
we selected those whose peak amplitude h matched that 
of the the observed cluster-template CCF, within a toler- 
ance Ah (Figure and used those data for our W v — a 
calibration. 

For CCF peak amplitudes comparable to those of our 
observations, we found that the relation between fitted 



CCF width W v and input a was always well-described 
by a linear relation whose parameters vary with template 
star and wavelength region. The value of the heliocentric 
velocity was not affected by low CCF peak amplitude, to 
within the errors. 

To obtain a realistic estimate of the errors in the ve- 
locity dispersion calibration, we defined an acceptable 
spread Ah about the measured CCF peak amplitude 
(^meas ^ A/i), and used only points lying in this inter- 
val. Here, h^ cas is the measured CCF peak amplitude for 
cluster C. The value of Ah defines how closely we require 
the peak amplitude of our simulated CCFs to match that 
of the CCF peak from cross-correlating the cluster data 
with the template stars. The uncertainty in the derived 
value of a is clearly related to the scatter in each fit (see 
Figure EJ, but it also depends on the value of Ah. For 
example, when Ah is large, the scatter 6 about the linear 
relation between W v and a is also large. 

For each combination of cluster, template star and 
spectral region, we adjusted the value of Ah until sw, 
the scatter in the CCF width, approached AW, the un- 
certainty we estimated from using f xcor to measure the 
same CCF width (e.g. Figure [BJ. In most cases, sw was 
somewhat larger than AW , even for arbitrarily small val- 
ues of Ah. Hence the uncertainty in the derived a values is 
not only due to uncertainty in measuring W v but includes 
a small additional component from the cross-correlation 
method itself. Hence, s a should give a realistic estimate 
of the uncertainty in each a measurement, and these are 
the values quoted in parentheses in Table [5] 

The derived velocity dispersions show some variation 
with template star, though where the CCF amplitude is 
highest (region Gr), the derived values of a are gratify- 
ingly similar. The F5 V template star HR 7749 gives sys- 
tematically lower values of a for both clusters, but only to- 
wards the blue, where we expect more contamination from 
early-type stars. For these reasons, and those outlined in 
Sect. 13.4.21 we regard the blue part of the G window as 
less reliable and prefer therefore the results from region 
Gr. The higher measured dispersion in Rr could be due 
in part to contamination from the broader lines seen in for 
example Sin in some early-type stars (uvespop database, 
Bagnulo et al. I2003J> . However, a velocity dispersion that 
varies as a function of wavelength is not necessarily un- 
physical but may arise naturally due to mass segregation, 
although we do not claim that this is the reason for the 
present discrepancy. 

As we argued in Sect. l3.4.H early- type supergiants are 
not likely contributors to the absorption spectra in the vi- 
sual region of either #23 or #34. We nevertheless carried 
out some tests with an A-type supergiant as template star, 
in order to further constrain the dependence of our results 
on the assumed template star. The A2 lab star HD 102878 
from the UVESPOP database (Bagnulo et al.l2"0"05) was used 
for this purpose. We found that while this star gave ac- 

6 We denote the RMS scatter in W v and a about the linear 
relation as s w and s a , respectively. 
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ceptable cross-correlation signals, indicating the same ra- 
dial velocities as for our standard set of templates, we were 
unable to use it for obtaining a robust a — W v calibration 
relation for our noisy data. Better results were obtained, 
however, for a hybrid template star which we constructed 
by adding HD 102878 and the K4 supergiant HR 6961 
in proportions 1:4. This hybrid star produced radial ve- 
locities and a fully consistent with those of our normal 
template stars. Hence, even if early-type stars could con- 
tribute to the continua, they would not contribute much 
to the CCF signal. 

3.4.4. Results for cluster #23 

In Tables and H3 we present the resulting velocity dis- 
persions and radial velocities for the different wavelength 
regions and template stars. 

For #23, we find a weighted average value of cr%z = 
32.5 ± 2.7 km s _1 in region Gr (Sect. I3.4.2|) . when using 
all 5 template stars, and (T 2 3 = 32.5±2.3 if we only use the 
two K supergiants (see discussion in Sect. 13.4.11 . The un- 
certainty here includes both the uncertainty in each mea- 
surement (s CT ; the values quoted in parentheses in Table 
03 see Sect. I3.4.3j) . and the smaller variation between re- 
sults for different template stars (i.e. the standard devia- 
tion of the a column in question). We use this method to 
estimate the total uncertainty in a in what follows. For 
region R we find 023 = 37.1 ± 3.4 km s _1 when using all 
stars and 023 = 36.7 ± 3.1 when using the K supergiants 
only. Combining the measurements from Gr and R we 
obtain <t 2 3 = 34.4 ± 3.8 km s _1 when using all stars, and 
023 = 34.1 ± 3.5 km s~ 1 when using only the K super- 
giants. 

Following the discussion in Sect 13.4.11 giving higher 
weight to the K super giants, and noting the better 
quality of the Gr data, we adopt as our best estimate 
(T23 = 32.5 ±2.5 km s _1 . We note that adopting a velocity 
dispersion of the magnitude implied from region R would 
produce a virial mass that is 25 ± 25% higher. 

We find heliocentric velocities «hciio,23 = 2859.7 
km s _1 in Gr and 2855.2 in R, with a scatter between 
the five template stars that in both cases is only ~ 1 
km s . Restricting to K supergiants only makes no signif- 
icant difference for Vhoiio,23- The output uncertainties on 
^heiio from f xcor are on the order of 4 km s _1 in Gr and 
8 km s -1 in R (see TableEJ). The small star-to-star scatter 
suggests that the uncertainty is dominated by the quality 
of our cluster spectra rather than by template star mis- 
matches. We note a systematic difference of ~ 4 km s _1 
between regions Gr and R, that however is within the 
uncertainties. We base our final value on the Gr data: 
Whdio,23 = 28 59.7 ±4.4. 

To investigate possible contamination from other clus- 
ters, we compared cross-correlation results for the 1.8- 
arcsec extraction with the narrower, 1-arcsec wide, extrac- 
tion around cluster #23 (see Sect. I2.2|> . The narrower ex- 
traction gives CCFs with poorer signal-to-noise ratio but 



very similar velocity dispersion in both Gr and R. Also, 
the use of hybrid templates made from HD 102878 (A 
supergiant) and HR 6961 (K supergiant) produced very 
similar results (023 = 30 km s _1 in both Gr and R). 

3.4.5. Results for cluster #34 

For cluster #34 (Tables |5] and EJ) the analysis was more 
straightforward. We measured 034 = 17.6 ± 2.0 km s _1 
in Gr and 17.4 ± 4.4, in R. Combining Gr and R yields 
CT34 = 17.5 ±2.8 km s _1 . 

For the radial velocities we find «hciio,34 = 2 7 94.6 ± 
2.2 km s _1 in Gr and 2791.7 ± 5.6 in R, i.e. we see a 
similar systematic difference between Gr and R as for 
cluster #23, but which again is within the uncertainties. 

Given the better quality of the Gr data, we adopt 
these results for 034 and i>hciio,34- 

3.5. Evidence for ISM gas flows from Na 1 lines 

We also carried out cross-correlations on the narrow spec- 
tral region 5874-5912 A, which is dominated by the Na I D 
doublet. For cluster #23 this was successful (see Fig. [SJ, 
but the wavelength region was too short to yield anything 
useful in the more noisy spectrum of #34. Nor could we 
obtain a reliable velocity dispersion calibration. The result 
would in any case be difficult to interpret, since we might 
also expect a contribution from ISM absorption. The he- 
liocentric velocity is well-determined and remarkably low. 
At 2845±7 km s _1 it is 25 km s" 1 lower than that inferred 
from region Gr along the same line of sight. Presumably, 
the Na 1 absorption we are seeing here is a combination 
of both stellar lines and cold interstellar medium moving 
at a speed of ~ 20 km s _1 along the line of sight to the 
cluster (cf. Schwartz and Martin 2004 ) . The emission lines 
also show evidence for a blue-shifted component (see Sect. 

3. 6. Cluster sizes 

We have measured the spatial extent of the two clusters 
using the HST images described and presented by Ostlin 
et al. (HHHJ). 

To measure the sizes of #23 and #34 we used curves of 
growth extending from to 5 Planetary Camera (PC) pix- 
els. One PC pixel has a size of 0".0455 which corresponds 
to 8.27 pc at the adopted distance of ESO338-IG04. We 
used observations in 4 passbands: F336W, F439W, F555W 
and F814W. 

The growth curves were compared to model PSFs gen- 
erated with Tiny Tim 7 (Krist and Hook 2003) that we 
convolved with Gaussian profiles and King profiles with a 
concentration parameter c = 2. 

Using the whole growth curve from to 5 pixels, the re- 
sulting sizes are well-constrained. We find the best-fitting 
effective half-light radii to be r e = 5.2 ± 1.0 and 5.6 ± 0.4 

7 http:/ /www. stsci.edu/software/tinytim 
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Table 5. Derived velocity dispersions in km s _1 for the two clusters, for each template star and spectral region. 
Errors are 1-sigma and are equal to the RMS scatter in the a — W v relation for each particular measurement, see text 
for details. 



ref. 


Type/ 




Spectral region (see 


Table 




star 


class 


G 


Gr 


Rb 


R 


Rr 


Cluster #23 


HR 6961 


K4lb 


33.4 (2.5) 


32.3 (2.2) 


27.4 (5.1) 


36.2 (3.1) 


38.6 (2.8) 


HR 7277 


Kli 


33.0 (3.2) 


32.8 (2.4) 


29.7 (6.7) 


37.2 (2.9) 


48.6 (7.8) 


HR 7706 


Kim 


32.5 (3.1) 


31.3 (3.1) 


29.3 (3.6) 


33.2 (4.8) 


41.1 (4.4) 


HR 7749 


F5v 


21.2 (2.7) 


32.0 (2.4) 


30.9 (3.5) 


39.7 (3.9) 


46.0 (6.2) 


HR 8299 


G5m 


35.4 (2.1) 


34.6 (3.3) 


39.5 (6.1) 


37.7 (2.4) 


43.8 (4.6) 


Hybrid 1 


K+A 


31.9 (3.4) 


29.9 (3.5) 


29.4 (5.0) 


30.2 (2.4) 


37.5 (13) 


Cluster #34 


HR 6961 


K4ib 


18.6 (2.4) 


17.4 (2.0) 


19.1 (2.8) 


15.8 (3.0) 


34.3 (19) 


HR 7277 


Kli 


18.2 (2.4) 


17.6 (1.7) 


16.2 (4.4) 


11.0 (5.6) 


18.2 (4.4) 


HR 7706 


Kim 


18.0 (2.5) 


16.8 (3.0) 


16.3 (3.5) 


15.8 (3.5) 


15.8 (8.0) 


HR 7749 


F5v 


9.1 (1.9) 


17.4 (2.0) 


16.0 (3.5) 


22.1 (3.3) 


55.8 (30) 


HR 8299 


G5m 


18.8 (1.7) 


18.1 (1.8) 


17.6 (8.2) 


18.1 (2.7) 


20.3 (32) 


Hybrid 1 


K+A 


18.5 (2.3) 


18.4 (2.4) 


15.6 (5.2) 


15.3 (2.4) 


18.5 (15) 



Note: (1) K4lb and A2Iab in proportion 4:1. 

Table 6. Derived heliocentric velocities in km s _1 for the two clusters, for each template star and spectral region. 
Errors are formal errors in Gaussian fitting to the CCF peak from f xcor. 



ref. 


Type/ 




Spectral region (see 


Table EJ 




star 


class 


G 


Gr 


Rb 


R 


Rr 


Cluster #23 


HR 6961 


K4lb 


2860 (7) 


2860 (4) 


2858 (7) 


2857 (7) 


2855 (7) 


HR 7277 


Kli 


2857 (7) 


2858 (4) 


2855 (6) 


2854 (7) 


2853 (7) 


HR 7706 


Kim 


2860 (8) 


2861 (5) 


2857 (9) 


2856 (9) 


2856 (9) 


HR 7749 


F5v 


2858 (8) 


2860 (5) 


2857 (8) 


2854 (9) 


2853 (8) 


HR 8299 


G5m 


2859 (7) 


2860 (4) 


2857 (7) 


2855 (7) 


2854 (7) 


Hybrid 1 


K+A 


2860 (7) 


2859 (4) 


2857 (8) 


2857 (7) 


2856 (7) 


Cluster #34 


HR 6961 


K4ib 


2794 (4) 


2795 (2) 


2793 (5) 


2792 (5) 


2789 (14) 


HR 7277 


Kli 


2792 (4) 


2793 (2) 


2791 (5) 


2790 (5) 


2789 (5) 


HR 7706 


Kim 


2794 (4) 


2796 (2) 


2793 (7) 


2793 (5) 


2793 (5) 


HR 7749 


F5v 


2793 (4) 


2794 (2) 


2793 (7) 


2790 (8) 


2785 (13) 


HR 8299 


G5m 


2794 (3) 


2795 (2) 


2793 (6) 


2793 (6) 


2792 (5) 


Hybrid 1 


K+A 


2794 (4) 


2795 (2) 


2792 (5) 


2792 (5) 


2791 (5) 



Note: (1) K4lb and A2Iab in proportion 4:1. 



pc for #23 and #34, respectively, when using a Gaussian 
profile. Since the difference in estimated effective radius is 
in general small when comparing Gaussians and different 
King model assumptions, r e is a quite model-independent, 
and hence robust, measure of the size of a cluster (Carlson 
and Holtzman lMfll Maraston et al. l2UU3jl . 

3. 7. Virial masses 

From the virial theorem, we can estimate the mass of a 
cluster if we know its velocity dispersion a and its size. 
Assuming the cluster is in virial equilibrium, its mass is 
given by 



M vir = ka 2 r m /G, 

(e.g. Spitzer ITW1 Smith and Gallagher I2TM1 McCrady 
et al. I2(JU3[) where r m is the half-mass radius. Taking 
k = 7.5, r m — |r e (Spitzer ri987fl . we find that the cluster 
masses are given by: 

M vh = 2.324x 103 ( 1 ^-t) 2 (^)m 
Specifically, for #23, 

M23 - L28Xl ° 7 U/kms-^ fe) M ° 
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and for #34, 

In Tabled we summarise the properties of clusters #23 
and #34, including the resulting dynamical mass-to-light 
ratios (M/Ly). The uncertainty intervals follows from a 
root-square addition of the uncertainties in a and r e . Boily 
et al. (J2Q05) discuss the evolution of r\ = k ■ r m /r e with 
time during the evolution of a cluster. For a cluster of the 
age of #23, i.e. ~ 5 Myr, Boily et al. PTO find r) » 10, 
which is what we have used above. For ages greater than 
50 Myr, however, they predict 77 « 20 for models with high 
initial density. 

3.8. Velocity structure in emission lines 

The slit covers a region of strong emission from ionised 
gas. In Figure El we show a spectral image centred on 
[O in] A5007 showing that the ionised gas motions along 
the slit are complex. The same features are seen in Ha. 
The wavelength of the peak of the Ha and [O ill] A5007 
lines vary by almost 100 km s _1 along the slit, but not in a 
monotonia fashion. Closer inspection reveals the presence 
of several velocity components at each position along the 
slit. For most positions at least three Gaussians appear to 
be needed to adequately reproduce the observed emission. 

Figure EH shows the Ha line profile at the location 
of #23. The line is very broad and clearly non-Gaussian, 
with peak intensity blue-shifted by 45 km s _1 with respect 
to the velocity inferred from the absorption lines in #23. 
We fit three Gaussians to the observed profile and find, 
one narrow blueshifted component (at —40 km s _1 with re- 
spect to the absorption line velocity), one redshifted com- 
ponent (at +20 km s _1 ), and one broad (a — 65 km s _1 ) 
component with nearly the same velocity as #23. We see 
the same pattern in both Ha and [O ill] 5007 A. 

Looking at the scan lines to the north-east of cluster 
#23 (above #23 in Fig. EJ, the line is split into two com- 
ponents: one blueshifted and one redshifted with respect 
to the centre of the cluster. For the scan lines immediately 
below cluster #23, the blueshifted component dominates. 

The blue- and red-shifted components suggest a bubble 
expanding with a velocity of ~40 km s _1 . If we take a 
characteristic radius of 1" (see Fig. [3] and Sect. 14. if) . the 
timescale for blowing such a bubble is ~5 Myr, which is of 
the same order as the age of #23. Our analysis here also 
shows that emission lines, being much broader than the 
cluster's absorption lines, cannot reliably be used to infer 
the dynamical mass of associated clusters (cf. Turner et al. 
2003J. A few scan lines further down (marked with ticks in 
Fig.[5J) we see high- velocity components, redshifted by 200 
km s _1 and 80 km s _1 . We finally note that cluster #34 
has a stellar absorption line velocity which is 65 km s _1 
less than the ionised gas from the same position on the 
slit. 



4. Discussion: Comparing cluster masses and ages 
from photometry and spectroscopy 

The ages and photometric masses of the cluster population 
in ESO 338-IG04 have been investigated by Ostlin et al. 
(2003). Here, we discuss the properties of cluster #23 and 
#34 in more detail, and compare with the virial mass 
estimates and age estimates based on our UVES spectra. 

For the sub-population of young (age < 20 Myr) clus- 
ters, Ostlin et al. H2003[) preferred a low metallicity for 
both the stellar and ionised gas components. This agrees 
with the nebular oxygen abundance of 12+log(0/H) = 8.0 
(Bcrgvall and Ostlin 12*0021) . Moreover, to reproduce the 
youngest clusters, the IMF needs to have a slope equal to 
or flatter than the Salpeter value(a = 2.35) and an upper 
mass limit M up > 60M Q . We shall refer to this preferred 
metal-poor model as the "standard" model 8 in the follow- 
ing sections. 

In the absence of dark matter, the photometric and 
dynamical masses should agree. If this is not the case, 
the reason may be sought either in the determination of 
the dynamical mass (because of mass segregation, or non- 
equilibrium), or the photometric mass (IMF, stellar pa- 
rameters). Below we discuss the photometric and dynam- 
ical masses of each cluster in this context. 

4.1. Cluster #23 

The standard model gives an age of 7 Myr, an internal 
reddening of Eb-v = 0-05 and a photometric mass of 
5 x 10 6 M Q , i.e. significantly lower than the virial mass. 
By adopting a steeper IMF slope (a = 2.85) or a low 
upper mass limit (M up = 20Af Q ), photometric masses in 
the range 2 to 6 x lO 7 A/ can be produced. 

The models of Zackrisson et al. (2001; Z01) that were 
used in Ostlin et al. I|2003[) include nebular emission. 
However, the nebular emission produced by each cluster 
is assumed to be included in the photometric aperture, 
which has a radius of 0".14, or 25 pc at the distance of 
ESO 338-IG04. This assumption appears to be valid for 
the majority of young clusters, but it does not hold for 
#23. This cluster alone accounts for just over 10% of the 
integrated starburst flux at 2000 A and has ionised a large 
bubble (radius ~ 1", i.e. nearly 200 pc), which, as we dis- 
cussed in Sect. 13.81 is expanding (Figs. 1^1 and mil see also 
Fig. 2 in Ostlin et al. 2003, or Fig. 7 in Hayes et al.EEEIl. 
Since ionised gas dominates the luminosity output when 
the cluster is young (Bergvall and Ostlin I2UU2|) . the total 
luminosity, and hence the photometric mass of #23 may 
be underestimated by at most a factor of about two (see 
Fig. 17 in Bergvall and Ostlin l2"()()2|l . For instantaneous 
bursts of age ~ 6 Myr, the correction to the mass is in- 
stead of the order +25%. The displaced nebular emission 
also means that the results of the photometric modelling 



Standard model: instant burst, IMF with slope a = 2.35 
(Salpeter) and mass range 0.08-120 Mq, gas covering factor of 
unity, metallicity of stars and gas: Z gas — 0.002, Z sta rs = 0.001. 
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Table 7 

for Galactic reddening according to Schlegel et al 



Properties of clusters 23 and 34. The colours, absolute magnitudes and mass-to-light ratios are corrected 

1998;), but not for internal reddening. 



Cluster no. 


#23 


#34 


Absolute magnitude M v 


-15.5 


-12.8 


Colour (t) — i) 


0.11 


0.86 


Effective radius r e (pc) 


5.2 ± 1.0 


5.6 ±0.4 


Heliocentric velocity Uheiio (km s _1 ) 


2859.7 ± 4.4 


2794.6 ± 2.2 


Velocity dispersion a (km s _1 ) 


32.5 ±2.5 


17.6 ±2 


Virial mass (1O 6 M ) 


13 ±3 


4.0 ± 1 


Mass-to-light ratio (M/Lv)q 


0.093 ± 0.02 


0.35 ±0.09 




#23 



2700 



2800 



2900 



3000 



Fig. 9. Two-dimensional spectrum around [O III] A5007 showing the complex velocity structure in the ionised gas. 
The x-axis is labelled with velocity in units of km s _1 . The y-axis is the location along the slit, and the tick marks 
have a spacing of one arcsecond. The position of #23 is indicated by an arrow, and two high velocity structures are 
marked by short dashes. 



for #23 are more uncertain than for most other young 
clusters in ESO 338-IG04. 

Indeed, recent Ha photometry with HST/ACS (Ostlin 
et al., in prep.), gives an integrated Ha luminosity for #23 
(within the same aperture as above) of ~ 10 32 W. For the 
electron density derived from our UVES spectra (Sect. 
13. 1[) . this implies a total ionised gas mass of ~ 10 4 M Q or 
about one thousandth of the dynamical mass. We caution 
that the density is derived from a larger region where the 
nebular component is dominated by the more extended 
gas around #23, but the value is sufficiently small that 
we can conclude that the mass of ionised gas in #23 is 
negligible compared to the dynamical and stellar masses. 
Moreover, the Ha emission line equivalent width for #23 
derived from the ACS images is only about ~ 50 A, much 
smaller than predicted for an age of ~ 7 Myr (Starburst99, 
Lcithcrcr et al. I1999fl . Hence, the ionising photons pro- 
duced in #23 leak out over a much larger volume than 
the size of the cluster. Much of the ISM originally present 
in this cluster is likely to have been removed already. 

Another reason for the discrepant masses could be that 
the cluster is not in virial equilibrium. This could be the 
case if the cluster is so young that it has not yet had time 
to virialise, or because supernova explosions have led to 



the prompt expulsion of gas from the cluster. In such cases 
its dynamical mass may be overestimated by a factor of 
up to 3, depending on the star formation efficiency (SFE) 
in the cluster, i.e. the mass fraction of gas initially present 
that was turned into stars (Bastian and Goodwin I2(JU6|I . 
In the case of #23 we would infer a SFE of ~ 40% for 
a stellar mass of 5 x 10 6 M Q . If we account for the dis- 
placed nebular emission, we obtain SFE > 50%, the exact 
amount depending on the dynamical state of the cluster 
prior to gas expulsion. If the fraction of gas lost is larger 
than 70% the cluster may dissolve completely, but even a 
cluster which remains bound will lose a few tenths of its 
initial stellar mass (Bastian and Goodwin 2006). Hence 
we expect #23 to remain bound but to lose a few tenths 
of its mass over the next ~50 Myr. 

Alternative constraints on the age and IMF can be 
obtained from the spectra presented in this paper. Our 
analysis of Balmer and He I absorption lines in Sect. 13.31 
indicated an age of 61 4 , Myr. The effect of adopting an 
IMF with relatively few massive stars is to increase the 
equivalent widths at young ages. A low upper mass limit 
(Mup > 30M Q ) may still be marginally consistent with 
the data, but an IMF slope significantly steeper than the 
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Fig. 10. Ha line profile for the two scan lines centred on #23, without background subtraction. The x-axis is labelled in 
velocity. The thin (black) solid line shows the observed spectrum. The (blue) dotted line shows the modelled spectrum 
which here is a sum of 3 Gaussian profiles, which are shown as (green) dashed, dash-dot and dash-dot-dot lines. The 
(red) dotted line at the bottom shows the residuals (observed — model spectrum). The vertical (black) long-dashed 
line shows the velocity of #23 determined from absorption lines through cross-correlation. This spectrum shows both 
that there is both a blueshifted and a redshifted component, and that there is a broad component with velocity close 
to the velocity of #23. 



Salpeter value predicts equivalent widths that are too high 
(see for example Fig. 11 in Gonzalez Dclgado et al. 1999). 

More evidence for the the presence of very massive 
stars comes from the detection of intermediate- width stel- 
lar Hen A4686 line emission with equivalent width 0.8 A. 
We see no obvious WR features but interpret this line 
as a signature of Of stars, which are very massive stars 
that may be on their way to becoming WR stars (Nota 
et al. 119961 Schaerer and Vacca I1998|l . This observation 
makes an IMF with a low upper mass limit improbable, 
but is consistent with a Salpeter IMF and an age < 5 Myr 
(Schaerer and Vacca 1 199*51 Leitherer et al. H999[) . 

We noted in Sect. 13.4.11 that the presence of strong 
calcium triplet absorption from #23 signals the presence 
of red supergiants. For a Salpeter IMF and low metallicity, 
the presence of RSGs indicates an age of > 5 Myr, whereas 
for an IMF with a low upper mass limit, the RSGs require 
that the cluster is older than 10 Myr. The simultaneous 
presence of RSGs and emission in He n A4686 is somewhat 



unexpected, but could be accounted for if the formation 
of the cluster extended over a few Myr, rather than being 
instantaneous. Similar scenarios have been presented for 
R136 (Massey and Hunter IT^gjl and NGC3603 (Pandey 
et al. l2000fl . Moreover, recent stellar models which include 
rotation predict longer WR lifetimes (Meynet and Maeder 
I2UU5|) . 

Fellhauer and Kroupa (2005) have proposed that 
the anomalously massive, intermediate-age cluster W3 in 
NGC 7252 may be the product of merging of several less 
massive clusters. Though the timescale for this to happen 
seems short in the case of #23, an earlier merger of young 
clusters could lead to the presence of stellar populations 
with different ages, and explain the simultaneous presence 
of Of stars and RSGs. 

It is also quite possible that stellar evolution in a clus- 
ter as massive as #23 is different compared to a lower- 
density environment. Numerical simulations of massive 
star clusters indicate that dynamical friction rapidly leads 
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to strong mass segregation and possibly merging of the 
most massive stars (e.g. Portegies Zwart et al. 1999). 
Massive-star mergers would likely have a sizeable impact 
on the early spectral evolution of such a star cluster. 

In conclusion, we find strong support for an IMF of 
=#=23 that extends to high masses (on the order of 100 
Mq) and has a slope not significantly steeper than the 
Salpeter value. The cluster's photometric and dynamical 
masses are in reasonable agreement, when allowing for the 
fact that most of its associated nebular emission lies out- 
side the photometric aperture, and the possibility that the 
cluster is not yet fully virialised. A more exotic possibil- 
ity would be that the cluster contains of the order of 50% 
dark matter. 

4.2. Cluster #34 

Cluster #34 is much older than both #23 and the current 
starburst. Our standard 8 metal-poor model (considered in 
Ostlin et al. 120031 based on the code by Zackrisson et al. 
2001; Z01), implies an age of ~ 1.5 Gyr, and photometric 
mass 2 x 10 7 M Q , almost an order of magnitude larger 
than the virial estimate ~4x 10 6 Mq. Compared to the 
photometry of #34, this model is discordant at the 1- 
sigma level. Varying the IMF slope a in the range 1.35 
to 2.85 gives comparable qualities of fit and a mass > 10 7 
Mq . A better fit can be obtained by adopting a model with 
significantly higher metallicity (Z = 0.008 to 0.04). This 
results in a lower age, 250 Myr (for Z = 0.04) to 500 Myr 
(Z = 0.008), and a somewhat lower mass of ~ 1O 7 M . For 
all assumed metallicities and IMF parameters, we find a 
rather high internal reddening, Eb-v ~ 0.25, which can 
be compared to the low values for the galaxy's starburst 
centre (E B _ V < 0.05; Ostlin et al. 2003) suggesting that 
#34 may lie on the far side of ESO 338-IG04. About one 
third of the other old and intermediate age clusters have 
similar implied reddenings. 

In addition to Z01, whose stellar evolutionary tracks 
come mainly from the Geneva group, we have used 
the spectral synthesis code pegase.2 (Fioc and Rocca- 
Volmerange 1997, 1999) which is mainly based on tracks 
from the Padova group. One effect of adopting a different 
set of tracks is that pegase.2 in general implies a lower 
reddening and higher age for #34. The best fit is again 
obtained for a rather metal-rich [Z — 0.008) model that 
gives an age of 1.2 Gyr, a reddening of Eb~v — 0.1 and 
mass 1.1 x 1O 7 M for a Salpeter IMF. However, pegase.2 
gives results that are consistent with the photometry for 
all metallicities in the range 0.002 < Z < 0.05, with re- 
sulting ages from 0.8 to 1.2 Gyr, and masses which differ 
by less than 20%. In summary, models based on a Salpeter 
IMF predict a mass of ~ 1O 7 M , irrespective of the set of 
stellar evolutionary tracks or metallicity used. 

All the mass estimates so far in this subsection assume 
a single power law IMF, including remnants but excluding 
the gas returned in the process of stellar evolution, which 
is assumed to have left the cluster. If we instead assume 



an IMF like that seen in the solar neighbourhood, with a 
flatter slope for stars with mass smaller than 0.5 — 1 Mq, 
the photometric masses quoted above are reduced by a 
factor of ~ 2. For instance, a Scalo98 (Scalo RM^i IMF 
yields ~ 5 x 10 6 Mq for the best fitting models. This is 
only slightly larger than the virial mass and such an IMF 
is also consistent with Galactic GCs (Chabrier and Mera 
ITS57ji . 

The work by Boily et al. 120051 suggests that for ages 
larger than 50 Myr, masses may be underestimated by 
a factor of two if the temporal evolution of r\ is not ac- 
counted for (see also Sect. 13 .7|) . However, this is only 
true for clusters with high initial mass surface density, 
and moreover, the evolution past 50 Myr is not known. 
Our #34 has a present mass surface density of ~ 2 x 10 4 
Mg/pc 2 (without any ^-evolution included), so strong ef- 
fects are not expected. 

Given the concordance of the photometric mass for 
different models and metallicities, we conclude that the 
dynamical mass of #34 is best explained with an IMF 
where the low-mass part is similar to that observed in the 
solar neighbourhood and in Galactic GCs. 

4.3. Implications for the globular cluster system in 
ESO338-IG04 

The two clusters studied in this paper are consistent 
with a normal IMF, extending to low stellar masses 
(M ~ O.IMq) and with a slope in the high mass regime 
(M > 1M Q ) close to the Salpeter value. Such IMFs will 
favour survival against cluster destruction mechanisms re- 
lated to stellar death and gas-expulsion. A few SSCs in 
other galaxies on the other hand seem to have rather odd 
IMFs, e.g. M82-F, for which a severe deficiency of low 
mass stars is implied (Smith and Gallagher 2001), making 
its long-term survival unlikely. 

All SSCs and GCs clusters line up nicely in a funda- 
mental plane (Walcher et al. ©05 ) suggesting that initial 
conditions regulate their properties (McLaughlin 2000). 
That young SSCs, including #23 and #34, follow the same 
trends suggests that the same basic formation mechanism 
applies both for globular clusters in ancient halos and for 
present-day mergers. Ostlin et al. (OZBR) used the clus- 
ter age distribution in ESO 338-IG04 to map its past star 
and cluster formation history, and found that star clusters 
seem to make up as much as several percent of the total 
stellar mass, in contrast to the value of 0.3% found for 
most galaxies (McLaughlin 2000). 

Our study demonstrates the massive nature of at least 
some of the star clusters in ESO 338-IG04. More velocity 
dispersion measurements in ESO 338-IG04 would be illu- 
minating. Though they are the brightest of the subpop- 
ulations of young and intermediate-age clusters, respec- 
tively, the objects studied in this paper are just two of the 
hundreds of clusters detected in this galaxy. Despite their 
impressive masses, #23 and #34 are apparently faint, 
with m v = 17 A and 20.1, respectively. Cluster #34 is the 
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faintest so far for which velocity dispersions have been 
derived from optical spectroscopy. With the exception of 
two very massive intermediate age clusters in NGC 7252 
(Bastian et al. I200b| . our clusters are also the most dis- 
tant yet probed. In this paper we have shown that velocity 
dispersion measurements of such faint clusters are indeed 
feasible, and we expect that clusters one magnitude fainter 
would be within reach with modest increase of integration 
time on VLT/UVES. Several more clusters in ES0 338- 
IG04 could be weighed by the same method. 

5. Conclusions 

We have presented high resolution spectra of two lumi- 
nous star clusters in the luminous blue compact galaxy 
ESO338-IG04 (Tololo 1924-416). The spectra have been 
cross-correlated with template stars observed using the 
same setup to determine each cluster's line-of-sight ve- 
locity dispersion. Using size estimates from Hubble Space 
Telescope images, we have used the velocity dispersions 
to determine the virial masses of the clusters. Our mass 
estimates have been compared to masses derived from 
HST photometry fitted to spectral evolutionary synthe- 
sis models. This comparison indicates that both clusters 
have rather normal IMFs, favouring their survival against 
internal disruption mechanisms. 

One of the clusters (#23) is young, with an age of 
^6 — 7 Myr and mass in excess of 10 7 M©, making it one 
of the most massive very young clusters known. We find 
evidence for the simultaneous presence of massive O-stars 
and red supergiants. The [O in] A5007 and Ha emission 
line profiles from the region surrounding cluster #23 indi- 
cate the presence of a bubble expanding at a velocity ~ 40 
km s _1 . In addition, we see signs of neutral gas flows along 
the line of sight from the cross-correlation analysis of the 
Na iD absorption lines. The IMF of this cluster shows no 
evidence for any deficiency of low-mass stars, but is con- 
sistent with a Salpeter slope over the whole mass range 
0.08-120 Mq. The inferred dynamical mass is larger than 
expected based on the luminosity of the cluster, which 
may indicate that the cluster is not yet virialised. 

The other cluster (#34) is older and has a mass of 
4 x 10 6 Mq. A combined analysis with the photometric 
data suggests that it has an age of between 0.3 and 1.4 
Gyr. The slope of the IMF at masses less than 1 M Q must 
be flatter than the Salpeter value for the virial and photo- 
metric masses to be consistent. Such an IMF is in agree- 
ment with that observed for old Galactic globular clusters. 
This work confirms that this cluster is a bona-fide young 
globular cluster. 

Acknowledgements. We thank Andrea Modigliani for useful 
suggestions on the reduction of UVES data. The authors ac- 
knowledge financial support from the Swedish research council 
and the Swedish National Space Board. We thank Matthew 
Hayes, Erik Zackrisson, Kambiz Fathi and Jan-Olov Persson 
for useful discussions. 



References 

Arp, H., Sandage, A., 1985, AJ, 90, 1163 

Bagnulo, S., Jehin, E., Ledoux, C, Cabanac, R., Melo, C, 

Gilmozzi, R., 2003, Messenger 114, 10 
Bastian, N., Saglia, R. P., Goudfrooij, P., Kissler-Patig, M., 

Maraston, C, Schweizer, F., Zoccali, M., 2006, A&A, 448, 

881, astro-ph/0511033 
Ba stian, N., Good win, S., 2006, MNRAS, 369, 9, 

astro-ph/0602465 
Bergvall N., Ostlin, G., 2002, A&A 390, 891 
Bertelli, G., Bressan, A., Chiosi, C, Fagotto, F., Nasi, E. 1994, 

A&AS, 106, 275 
Bica, E., Alloin, D., 1986, A&A 162, 21 
Bruzual, G., Chariot, S., 2003, MNRAS 344, 1000 
Boily, CM., Lancon, A., Dieters, S., Heggie, D.C., 2005, ApJ, 

620, L27 

Carlson, M. N., Holtzman, J. A., 2001, PASP, 113, 1522 
Chabrier, G., Mera, D., 1997, A&A 328, 83 
Diaz, A. I., Terlevich, E., Terlevich, R. 1989, MNRAS, 239, 325 
Drilling, J. S., Landolt, A. U. 2000, in Cox A. N., ed., Allen's 

Astrophysical Quantities, 4th edn., Am. Inst. Phys., New 

York 

Drinkwater, M.J., Gregg, M.D., Hilker, M., 2003, Nature, 423, 
519 

Dolphin, A.E., 2000, PASP 112, 1397 

Fellhauer, M., Kroupa, P., 2005, MNRAS, 359, 223 

Fioc, M., Rocca-Volmerange, B., 1997, Afc A 326, 950 

Fioc, M., Rocca-Volmerange, B., 1999, |astro-ph/9912179| 
(pegase.2) 

Gilbert, A. M., Graham, J. R. 2002, in IAU Symp. 207, 
Extragalactic Star Clusters, ed. E. K. Grebel, D. Geisler 
& D. Minniti 

Gonzalez Delgado, R. M., Leitherer, C, Heckman, T., Cervino, 

M. 1997, ApJ, 483, 705 
Gonzalez Delgado, R. M., Leitherer, C. 1999, ApJS, 125, 479 

(GDL99) 

Gonzalez Delgado, R. M., Leitherer, C, Heckman, T., 1999, 
ApJS, 125, 489 

Gray, D. F. 1981, ApJ, 251, 152 

Gray, D. F. 1984, ApJ, 281, 719 

Gray, D. F., Toner, C. G. 1987, ApJ, 322, 360 

Hayes, M., Ostlin, G., Mas-Hesse, J.M., Kunth, D., Leitherer, 
C, Petrosian, A., 2005, A&A 438, 71, astro-ph/0503320 

Ho, L. C, Filippenko, A. V., 1996a, ApJ, 466, L83 

Ho, L. C, Filippenko, A. V., 1996b, ApJ, 472, 600 

Kharchenko, N.V., Piskunov, A.E., Scholz, R.-D., 2004, 
Astron. Nachr., 325, 439 

Krist, J., Hook, R., 2003 "The Tiny Time User's Guide", ver- 
sion 6.1, |http://www.stsci.edu/software/tinytim| 

Larsen, S. S., Brodie, J. P., Elmegreen B. G., Efremov, Y. N., 
Hodge P. W., Richtler, T. 2001, ApJ, 556, 801 

Larsen, S. S., Brodie, J. P., Hunter, D. A. 2004, ApJ, 128, 2295 

Leitherer, C, Schaerer, D., Goldader, J.D., et al., 1999 ApJS, 
123, 3 

Maoz, D., Ho, L. C, Sternberg, A. 2001, ApJ, 554, L139 
Maraston, C, Bastian, N, Saglia R.P., Kissler-Patig M., 

Schweizer, F., Goudfrooij, P., 2004, A&A 416, 467 
Massey, P., Hunter, D.A., 1998, ApJ, 493, 180 
McCrady, N., Gilbert, A. M., Graham, J. R., 2003, ApJ, 596, 

240 

McLaughlin, D.E., 2000, ApJ 539, 618 

Mengel, S., Lehnert, M. D., Thatte, N., Genzel, R., 2002, A&A, 
383, 137 



Ostlin et al.: Dynamical masses of globular clusters in ESO338-IG04 



17 



Meurer, G.R., Heckman T.M., Leitherer, C, Kinney, A., 

Robert, C, Garnett, D.R., 1995, AJ, 110, 2665 
Meynet, G., Maeder, A., 2005, A&A, 429, 581 
Nidever, D. L., Marcy, G. W., Butler, R. P., Fischer, D. A., 

Vogt, S. S., 2002, ApJS, 141, 503 
Ochsenbein, F., 1980, Bull. Inf. CDS, 19, 74 
Ochsenbein F., Bauer P., Marcout, J., 2000, A&AS, 143, 221 
Nota, A., Pasquali, A., Drissen, L., Leitherer, C, Robert, C, 

Moffat, A. F. J., Schmutz, W. 1996, ApJS, 102, 383 
Origlia, L., Leitherer, C, Aloisi, A., Greggio, L., Tosi, M. 2001, 

AJ, 122, 815 

Ostlin, G., Bergvall, N., Ronnback, J., 1998, A&A, 335, 85 
Ostlin G., Amram P., Bergvall N., Masegosa J., Boulesteix J., 

Marquez I., 2001, A&A, 374, 800 
Ostlin, G., Zackrisson, E., Bergvall, N., Ronnback, J., 2003, 

A&A, 408, 887, astro-ph/0306522 
Pandey, A.K., Ogura, K., Sekiguchi, K., 2000, PASJ 52, 847 
Portegies Zwart S.F., Makino, J., McMillan, S.L.W., Hut, P., 

1999, A&A 348, 117 
Salpeter, E., 1955, ApJ, 121, 161 

Scalo, J. 1998, in ASP Conf. Ser., 142, The Stellar Initial Mass 

Function, 201 
Schaerer, D., Vacca, W. D., 1998, ApJ, 497, 618 
Schlegel, D.J., Finkbeiner, D.P., Davis M., 1998, ApJ, 500, 525 
Schwartz, C. M., Martin, C. L., 2004, ApJ, 610, 201 
Simkin, S.M., 1974, A&A, 31, 129 

Smith, L. J., & Gallagher, J. S., III. 2001, MNRAS, 326, 1027 
Spitzer, L., Jr., 1987, Dynamical Evolution of Globular 

Clusters, Princeton Univ. Press, Princeton, NJ 
Tonry, J., & Davis, M., 1979, AJ, 84, 1511 
Turner, J. L., Beck, S. C, Crosthwaite, L. P., Larkin, J. E., 

McLean, I. S., Meier, D. S., 2003, Nature 423, 621 
Underhill, A. 1995, ApJS, 100, 461 
Vacca, W. C, & Conti, P. S., 1992, ApJ, 401, 543 
Walcher, C.J., van der Marel, R.P., McLaughlin, D., et al., 

2005, ApJ 618, 237 
Werk, J. K., Jangren, A., Salzer, J. I., 2004, ApJ, 617, 1004 
Whitmore, B.C., 2003, in "A Decade of Hubble Space 

Telescope Science", Eds. Livio M., Noll K., Stiavelli M., 

Space Telescope Science Institute symposium series, Vol. 

14, p. 153, Cambridge University Press, astro-ph/0012546 
Whitmore, B., Heyer, I., Casertano, S., 1999, PASP, 111, 1559 
Whitmore, B. C, Zhang, Q., Leitherer, C, Fall, M. 

S.,Schweizer, F., Miller, B. W., 1999, AJ, 118, 1551 
Worthey, G., Faber, S. M., Gonzalez, J. J., Burstein, D. 1994, 

ApJ, 94, 687 

Zackrisson E., Bergvall N., Olofsson K., Siebert A., 2001, A&A 
375, 814 (Z01) 



